Evolution of Axionic Strings and Spectrum of Axions Radiated from Them 



Masahide Yamaguchi 1 , M. Kawasaki 2 *, and Jun'ichi Yokoyama 3 ^ 
1 Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 
2 Institute for Cosmic Ray Research, University of Tokyo, Tanashi 188-8502, Japan 
3 Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

(February 1, 2008) 



OS 
OS 

c 

m 

(N 
> 



oo ■ 

Oh! 
CD . 

43 ■ 



X 



Cosmological evolution of axionic strings is investigated by numerically solving field equations of 
a complex scalar field in 3+1 dimensions. It is shown that the global strings relax to the scaling 
solution with a significantly smaller number density than the case of local strings. The power 
spectrum of axions radiated from them is calculated from the simulation data, which is found to 
be highly peaked around the Hubble scale, and a more accurate constraint on the Peccei-Quinn 
breaking scale is obtained. 



PACS: 98.80.Cq, 11.27.+d, 14.80.Mz 

The global Peccei-Quinn(PQ) J7(l)-symmetry is in- 
troduced into the standard model of particle physics in 
order to solve the strong CP problem of the quantum 
chromodynamics (QCD) JlJ]. The axion appears as a 
Nambu-Goldstone boson in consequence of the sponta- 
neous breakdown of such a global symmetry, and acquires 
a mass through the instanton effect at the QCD scale [||. 

This breaking scale, f a , is constrained by the terrestrial 
and astrophysical experiments as well as by cosmological 
considerations. The most stringent lower bound has been 
obtained from SN 1987A as /„ > 10 10 GcV @. On the 
other hand, the upper bound is given by the condition 
that the present energy density due to the axion should 
not overdose the universe. 

In the very early universe, the i7(l)pQ-symmetry is ex- 
pected to be restored due to the high temperature effects. 
As the universe cools down, it spontaneously breaks down 
and global strings (axionic strings) are formed H . Subse- 
quent evolution of global strings has been less examined 
than that of local strings, and the result for the latter has 
been applied to axionic strings without direct numerical 
verification. In particular, global strings have been as- 
sumed to obey a scaling solution like local ones, in which 
the energy density of infinite strings is given by 



P = ^/t 2 



(1) 



where \i is the average energy per unit length of strings 
and £ is a constant. It has been found that £ ~ 13 for 
local strings in the radiation dominant era || . 

But there is a crucial difference between global and lo- 
cal strings; that is, the former strongly couples with the 
associated Nambu-Goldstone fields. As a result, global 
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strings have the following prominent features: (a) These 
fields carry most of the energy rather than string cores, 
(b) long-range forces proportional to the inverse separa- 
tion work between strings, and (c) the dominant energy- 
loss mechanism of loops is radiation of Nambu-Goldstone 
bosons || . Therefore it is not a trivial problem whether 
the global string network obeys the scaling solution like 
local ones. In fact, a deviation from scaling property is 
observed in 2+1 dimensions [Mpl, where "strings" cannot 
intercommute but only pair-annihilate. Of course, since 
strings intercommute with the probability of order unity 
in 3+1 dimensions, the results in 2+1 dimensions do 
not apply directly. In investigating the property of a 
global string, the effective action is often used, which is 
valid when massive modes are sufficiently damped. For a 
global string we have to use the Kalb-Ramond action [^0) , 
which is much harder to deal with than the Nambu-Goto 
action for a local string. Hence in the present Letter, 
instead of solving equations of motion derived from the 
Kalb-Ramond action, we numerically solve the evolution 
equation of the complex scalar field to trace the evolution 
of the global string network and calculate the spectrum 
of axions radiated from them. (If the reheating tempera- 
ture after inflation is lower than the symmetry breaking 
scale, global strings and radiated axions are washed out 
and relic axions come only from coherent oscillations of 
the zero mode Jll]]. Here we assume that this is not the 
case.) 

Axions are massless when they are emitted from strings 
and hence they behave like radiation. Later around the 
QCD scale, they acquire a small mass, m a , through QCD 
instanton effects. Since at present kinetic energy of ax- 
ions is much smaller than the rest mass, the present en- 
ergy density of axions, p a , is given by p a — m a n a with n a 
being their number density. Thus, in order to estimate 
the present energy density of axions, we must transform 
the energy density of radiated axions into the number 
density by use of the spectrum. However, the spectral 
shape is still in dispute. Davis, Shellard, and co-workers 
insist that the spectrum has a sharp peak at the Hubble 
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horizon scale ]12|-|14|. On the other hand, Sikivie, and 
co-workers claim that it is proportional to the inverse 
momentum . We clarify the shape of this spectrum 
in realistic situations where global strings evolve accord- 
ing to the scaling solution; that is, we identify kinetic 
energy of axions emitted from strings and elucidate the 
spectrum by Fourier transforming them. Hence our ap- 
proach is free from the ambiguities associated with the 
ad hoc choice of the initial configuration of strings. 

We consider the following Lagrangian density for a 
complex scalar field &(x), 



V cS [$] 



with the potential V c « [<&] given by 



(2) 



(3) 



Hereafter we set A = 1.0 for brevity. For T > T c = y/3r), 
the potential V e g has a minimum at $ = 0, and the U(l) 
symmetry is restored. On the other hand, new minima 
I $ | min = lyl ^ (T/T c ) 2 appear and the symmetry is 
broken for T < T c . In this case the phase transition is of 
second order. 

In the flat Friedmann universe, the equation of motion 
is given by 



$(s) + 3H$(x) 



1 



a(t) 



rVH(x) = ~V; S {<P] , (4) 



where the prime represents the derivative d/d& and a(i) 
is the scale factor. The Hubble parameter H = a(t)/a(t) 
and the cosmic time t are given by 



H 2 = 



Sn 



3M|j 30 



9*T 4 



t = 



1 
2H 



(5) 



where Mpi is the Planck mass, is the total effective 
number of relativistic degrees of freedom, and radiation 
domination is assumed. We define the dimensionless pa- 
rameter £ as 



< 



m(. 



\ 167T 3 <7*77 2 



1/2 



(6) 



In our simulation, we take £ = 10 and g„ = 1000, which 
corresponds to rj ~ 10 16 GeV, but the essential result is 
independent of this choice. 

We take the initial time U = t c /A and the final time 
tf = 75 ti = 18.75 t c , where t c is the epoch T = T c . Since 
the U(l) symmetry is restored at the initial time t = U, 
we adopt as the initial condition a thermal equilibrium 
state with $'s mass equal to the inverse curvature of the 
potential at the origin. The scale factor a(t) is normalized 
as a(U) = 1. 

We perform the simulations in 256 3 lattices with the 
physical lattice spacing fe p h ys = v3tja(t)/25. The time 



step is taken as St = O.Olt^ . Hereafter the subscript 
"phys" is omitted. The box size is nearly equal to the 
Hubble horizon H^ 1 and a typical width d ~ 1.0/(y/2rj) 
of a string is twice as large as the lattice spacing at the 
final time tf. We simulate the system from 10 different 
thermal initial conditions using the second order leap- 
frog method and the Crank-Nicholson scheme. We use 
the same method to identify a string core as in our pre- 
vious work H . More details will be published elsewhere 
]l6[ . The zero temperature potential V^ff[$,T = 0] is 
used after t — 20ti so that the axion is identified as in 
Eq. (||) below. We find that after some relaxation period 
the energy density of strings p is given by 



P = ^/t 2 



(J) 



where £ ~ (1.00 ± 0.08) irrespective of time and fi = 
2-Kiri 2 \a.{t/ (d^ 1 / 2 )) is the average energy per unit length 
of strings. Therefore, we can conclude that global strings 
network relax into scaling regime. In Fig. |l|, we show time 
development of £. 
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FIG. 1. Symbols(D) represent time development of £. The 
vertical lines denote a standard deviation over ten different 
initial conditions. 

Now we turn to the spectrum of axions emitted from 
axionic strings under the situation where they follow the 
scaling solution. If we represent the complex field <I>(i, x) 
in terms of the radial mode %(t, x) and the axion field 
a(t,x) as 



Ht,x) = 



V - 



X(t,x) 
V2 



exp 



ia(t, x) 



(8) 



the kinetic energy density of axions is given by 



1 



-d(t, x) 2 = , 
2 V ' ' |$0,x)| 4 



2 



-Im$(i, x)Re$(t, x) + Re$(i, x)lm.3>(t, x) 



(9) 



Since emitted axions damp like radiation, the average 
energy density of axions radiated in the period between 
ti and t 2 , p[ti,t2], is given by 

p[tx,ti] = y J d 3 xp[ti,t 2 ;x] 

= ^Jd 3 x ^a(t 2 ,x) 2 -^a(h,x) 2 (^- 



d 3 k 

V J (2tt) 3 

where V is the simulation volume and a(t, x) = 

I (f^MfcW exp(ifc-a;). 

To avoid contamination of string cores to the spectrum 
of emitted axions, we divide the simulation box into 8 
cells and stock the field data of a cell if there are no 
string cores in that cell between t\ and t%. Over all such 
cells, we average power spectra of kinetic energy of axions 
obtained through Fourier transformation. We follow the 
above procedure between t\ = 65ti and t 2 = 75t{, whose 
result is depicted in Fig. |[ 

Spectrum of axions emitted between L = 65L j and t=75l 
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FIG. 2. Filled dots represent the power spectrum, pk, 
which is already averaged over the direction of k and multi- 
plied by the phase-space factor as defined in Eq. (jj^) . It is ob- 
tained from cells with no string cores. The dotted line denotes 
its standard deviation. Open dots represent the spectrum ob- 
tained by averaging over all cells including string cores. Bins 
are cut every 55k. k = 645k corresponds to string cores. 



As is seen there, the spectrum, pk, is peaked around 
the horizon scale and decays exponentially for higher mo- 
menta. Note that it is contributed from both infinite 
and loop strings. For comparison, we have also given the 
power spectrum averaged over all cells, some of which 
include string cores. Note that effects of string cores lift 
the spectrum at high momenta and make it flatter. 

Defining the average inverse momentum by 



f°o dk pfc 
Jp 27T 7 k 

Jo 



dk „ 
2^Pk 



(11) 



we find k- 1 ~ (0.25 ± 0.18) x (i/27r). 

Finally, we want to estimate energy density of relic 
axions radiated from axionic strings. Though dynamic 
range of our simulation is limited to h\{t/d) ~ 5, we 
extrapolate our results up to the cosmological scale with 
ln(t/d) ~ 75. From the scaling property, we can estimate 
the energy density of axions radiated from strings from 
t = t* to r = t, 



Pa(t) 



2 -4 



i 



In 



dr 



(12) 



where t* is the cosmic time when axionic strings began 
radiating axions and we denote the breaking scale, 77, 
by f a . Multiplying it by the average inverse momentum 
/c _1 (t) = r/(27re), the number density of relic axions is 
given by 



n a (t) = 2nf^ 



1 



In 



de /2 



k- l (r)dT 



t e 



hi 



de /2 



(13) 



As the temperature cools down to the QCD scale and 
the axion acquires a non- vanishing mass, a network of 
domain-walls bounded by strings is created and walls 
start to dominate the dynamics of the system at t = t w 
given by p), 



1.7 x 10~ 6 A 2 



6 x 10" b eV 



fa 



10 12 GeV 



Ac 



QCD 

60 



0.5 



sec . 



-1.64 



(14) 



where A is a constant of order unity which describes un- 
certainties at the QCD phase transition, 



A = 10 



±0.5 



( Aqcd A 
V200MeVy 



6 x 10" 6 eV 

-0.65 , M 



fa 



10 12 GeV 

0.41 



0.82 



QCD 

60 



(15) 



with m a being the axion mass at zero temperature, Aqcd 
the energy scale of the QCD phase transition, and Aqcd 
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an effective number of massless degrees of freedom at 
that time. These domain-walls bounded by strings also 
decay by emitting axions if the QCD anomaly factor is 
unity fTj| . (Otherwise, these domain walls would rapidly 
over-dominate the universe.) 

As is seen from Eq.(|l3|), the dominant contribution to 
the present axion density is those radiated just before 
wall domination. The present density parameter, f2 s , of 
relic axions due to emission from strings is given by 



2.7 



h~ 2 A 



la 



10 12 GeV 



1.18 



(16) 



where h is the Hubble constant in units of 
100 km/sec/Mpc. 

With £ - (1.00 ± 0.08) and e" 1 - (0.25 ± 0.18), VL S is 
given by 

/ f \ 118 
n.~(0.68±0.46)fc- a A^ i5i ^J . (17) 

The condition that f2 s < 1.0 constrains the breaking scale 
of PQ-symmetry f a as 



f a < (1.39 ± 0.79) h 1 - 1 x 10 12 GeV , 



(18) 



which reads f a < (0.19-1.5) xlO 12 GeV for h = 0.5-0.8. 
Our constraint lies in between those found in [fl3f and in 
iPI , because the spectral shape agrees with the former 
but £ = 1 is substantially smaller than the case of local 
strings, £ = O(10), which has been adopted in |l3| . 

There are two additional contributions to the present 
energy density of axions besides that from strings con- 
sidered above. One is that from decay of a domain-wall 
connected by strings |18] and estimated as 



7/1 



fa 



10 12 GeV 



1.18 



(19) 



from which we obtain f a < 7 85 h 17 x 10 12 GeV with 
7 a factor of the order unity. Note that 7 also becomes 
smaller of the order 10 than the previous estimation. The 
other comes from the coherent oscillation of the axion 
field iflQfl and is estimated as 



0. = 1.5 x 10 



±0.4 



fa 



10 12 GeV 



1.18 



(20) 



which leads to f a < (0.33 - 1.5) x 10 12 GeV. These con- 
tributions may be comparable with or larger than axions 
radiated from axionic strings. Taking all these into ac- 
count, it is safe to say that f a < 10 11-12 GeV. 

In summary we have solved equations of motion of 
a complex scalar field to clarify cosmological evolution 
of axionic strings. As a result we have confirmed these 
global strings relax into the scaling solution but with a 
number density much smaller than local strings would 



occupy. We have also calculated the energy spectrum of 
axions generated from these axionic strings, which agreed 
with the result of Davis and Shellard |p|| . The constraint 
on f a turned out to be in between those obtained in [T^ ] 
and in p5[ . 
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